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In this paper, we investigate one-loop radiative corrections to the Casimir energy in the presence 
of two perfectly conducting parallel plates for QED theory within the renormalized perturbation 
theory. In fact, there are three contributions for radiative correction to the Casimir energy, up to 
order a. Only the main term, which is two- loop, has been computed by Bordag et. al (1985). Here, 
we consider corrections due to two remaining terms, i.e., photonic and fermionic one- loop correction 
resulting from renormalized QED Lagrangian. Our results show that only the fermionic loop has a 
minor correction. 



I. INTRODUCTION 

£NJ ■ The Casimir effect is the physical manifestation of the change in the zero point energy of a quantum field for 
different configurations. The zero point configuration refers to one in which there does not exist any on-shcll physical 
excitation of the field. 

In 1948 Casimir predicted the existence of this effect as an attractive force between two infinite parallel uncharged 
perfectly conducting plates in vacuum [l|. This effect was subsequently observed experimentally by Sparnaay in 
1958 [H (for a general review on the Casimir effect, see Refs.[5, jjf). Similar measurements have been done for other 
geometries, and their precisions have been greatly improved [5j-|l0(. The manifestations of the Casimir effect have 
been studied in many different areas of physics. For example, the magnitude of the cosmological constant has been 
estimated using the Casimir effect fllTll3{ . This effect has been also studied within the context of string th eory 
[T3 |. It has been investigated in connection with the properties of the space-time with extra dimensions fl5l - ll7j ]. 
The majority of the investigations related to the Casimir effect concerns with the calculation of this energy or the 
consequence forces for different fields in different geometries, such as parallel plates [l|, [Hj], cubes (l9l - l27j . cylinders 

Hi-Ill, and s P herical geometries [H, UMl • 

Although the Casimir effect has been known for nearly 60 years, the question of what are the leading radiative 
corrections to this effect is still surprisingly a subject of discussion. The first endeavors to compute the radiative 
corrections to the Casimir energy were reported in a paper by Bordag, Robaschik, and Wieczorek (BRW) [34 1 . T here 
exist many works on the radiative corrections to the Casimir energy for various cases (see for example [34j-|44|). In 
the case of a real massive scalar field, Next to Leading Order (NLO) correction to the energy has been computed in 
27 . l40l - 449l | . Moreover, the two- loop radiative corrections for some effective field theories have been investigated 
3a-l38]]. Bordag and his collaborators have calculated radiative correction to the Casimir energy due to one of 



the three related terms of order of a, ^/\/\/^ , in the presence of two perfectly conducting parallel plates for QED 

theory. In this view point, the photon propagator satisfies boundary conditions on the plates, while the plates are 
not transparent to the electrons. They found the correction Eq ' = 25 g ° m 4 to the popular leading term of Casimir 

energy = — 7 * 0a3 , where a is the distance between plates and m is the electron mass. In 1998 this result with 
another approach has been reported (35| . In the framework of renormalized perturbation theory for QED there are 
three vacuum bubbles of order of a. Up to now, all the papers on the Casimir effect that we are aware of, have 

not been calculated two of these diagrams to the radiative corrections, namely: photonic loop ^ resulting from 



Electromagnetic field and fermionic loop x from spinor field. 

The primary purpose of this paper is to directly calculate radiative correction to the Casimir energy resulting 
from one-loop corrections namely: one-loop photon and one-loop fermion, in the framework of the renormalized 
perturbation theory for QED theory. These corrections arc of order a. In order to do this, we use Green's functions 
in the presence of plates for Electromagnetic field with Dirichlct boundary condition and for spinor field with MIT 
bag boundary condition as propagators. Our main regularization is dimensional regularization. However, during the 
calculations we need cut-off regularization to eliminate infinities. 
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Our approach in the calculation of radiative corrections to Casimir energy is the most direct one. In this way we 
subtract two infinite energies: one relate to presence and the other without presence of two plates. We justify both 
of them to such a way our physical result is obtained. 

We organized the paper as follows: In Section [TT] we review briefly renormalization of Quantum Electrodynamics. 
In Section HTT1 using analogies between an Electromagnetic field and a massless scalar field, photonic loop correction is 
considered. We use the Dirichlet boundary condition on the two plates. In Section HVl we directly calculate radiative 
correction to the Casimir energy resulting from fermionic loop where MIT bag boundary condition as constraints on 
both of the plates is considered. In Section 5 we summarize our results and state our conclusions. 

II. RENORMALIZATION OF QUANTUM ELECTRODYNAMICS 

In this section we briefly review systematics of renormalization for QED theory (see for complete details 50]). The 
original QED Lagrangian is 

£qed - -\(F^f +i>(ip-m )i>-e^^At l . (1) 

i i_ 
By replacing ip = z 2 2 ip r and = z 3 2 A!f , it becomes 

£qed = -^z 3 (F^) 2 + z 2 ipr(i@ - m )i/j r - e z 2 z 3 2 tp r ^ip r A!^ , (2) 

where e is the physical electric charge and z 2 and z 3 are the field-strength renormalizations for ip and A M respectively. 
We define a scaling factor z\ as follows: 

i 

ezi = e z 2 z 3 2 . (3) 
We can split each term of the Lagrangian into two pieces as follows: 

£qed = - \{Fl v ) 2 + M*fi-™)A - ei; rl ^ r A^ - ^S 3 (FD 2 + M^~S m )A - eSifapilfrA?, (4) 

1_ 

where £3 = z% — 1, 62 = ^2 ~~ 1, $m = 2:2^0 ~~ m and Si = Z\ — 1 = (^-)z2Z^ — 1 are counterterms and m is the physical 
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mass of the electron. Now, the Feynman rules for this Lagrangian are 




(5) 




(6) 



~k 



— — (V v -(i-0— 

k 2 + ie\ k 2 



(7) 



(8) 



p — m + ie 



(9) 



i(p5 2 - 5 m ). 



(10) 



Each of the four counterterms must be fixed by renormalization conditions. For QED theory these conditions are (see 
for instance [50} ) 



»-^sj(n^sss, v = *n^( g )^( 3 ^fc 2 -fc^)n(fc 2 ) 



(ii) 




-»E(p) 



(12) 
(13) 



amputated 



In the above equations ill(fc 2 ) defines the sum of all 1PI (1-particle-irreducible) insertions into the photon propagator. 
Up to leading order in a it becomes 



n ( fc2 ) = 7~rT2" I dx E(L_a) —SasCl - x). 



(4tt)§ Jo f ^ 2 - 
m 2 — x(l — x)k 2 



(14) 
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Moreover, —iY,(p) denotes the sum of all 1PI diagrams with two external fermion lines. To leading order in a it is 



2™ /•! 



em 



(47r)2 J 



r(2- 1 



iS(p) = - l _ T / da;- i ^ r5 - 7 [(4-e)m-(2-c)a 3 »], (15) 



(1 — x) 2 m 2 + xfi 2 — x(l — x)p 



with e = A — d. Also, T fl (p',p) denotes the sum of vertex diagrams, more accurately 

FV-P) = 7^i(fe 2 ) + ^^F 2 (fc 2 ), (16) 

with iq and are unknown functions of k 2 called form factors and = f [7^,7"]- To lowest order, iq = 1 and 
F 2 = 0, we have T^ 1 = 7**. By using Eqs. ([T3 |I .([TB |1 and (fTo) . up to leading order in a, the counterterms are derived 
as follows: 

-p 2 r 1 r(2 - -) 

S z = ~ ~ I dx 7 ± 1 08x(l-x), (17) 



(47r)2 Jo (m 2 ) 



—em 



2™ /•! 



(47r) 2 Jo 



r(2-| 



*™ = 77-T / " . 2 _* ( 4 - 2a; -a;)), (18) 



(1 — x) 2 m 2 + x/i 



s * - f di 7 — r(2 ' l) v -« (< 2 - * - l ( i^T> - 21 - E<1 - *» I • (19 > 

(1 — x) 2 m? + xp? 



h . ^1 /V-.)f ^ (20) 

(47r)#y V((l-z) 2 m 2 + ^ 2 ) 2 -^ 2 

r(3 '^ } ^[2(l-4z + z 2 )- £ (l-,)> 2 Y (21) 



((1 — z) 2 m 2 + zfi 2 ) 3 ' 

where /1 is a photon mass to control infrared divergences which will finally go to the zero. According to the above 
discussion three vacuum bubbles contribute to the Casimir energy: j>, ), ^\/\^ . Two first diagrams arise 

from Eqs. (J5J) and (|10p. Bordag et. al have computed the last one. In the next sections we will consider the other 
two. 



III. PHOTONIC LOOP 



In this section, we calculate NLO radiative correction to the Casimir energy due to the photonic loop. We use 
Dirichlct boundary condition on the two parallel perfectly conducting plates in (3+1) dimensions. Obviously in the 
presence of the two plates, propagators automatically incorporate the boundary conditions and are position dependent. 
The contribution of one-loop photon to the vacuum energy in the interval [-^p, §] is 

/a/2 /" a /2 
d 3 x(n\H I \n) = i I d 3 x + 0(a 2 ), (22) 

-a/2 J-a/2 

using Eq.© it becomes 

/a/2 
d 3 x G B {x,x)[-i(g^k 2 - kn v )5z], (23) 
-a/2 

where Gb(x,x ) is the propagator of Electromagnetic field in the bounded space. For overall consistency, we use 
dimensional regularization to control ultraviolet divergences, and a photon mass fj, to control infrared divergences. 
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Using analogies between an Electromagnetic field and a massless scalar field, photon propagator is considered as 



G B (x,x ) 



,l*-*k± e- iw ^- t 'h- ik ^-( x± - x ' ± ^ sin(fc„(z + §)) sin(fc n (z' + §)) 



2tt J (2?r 



id-2 



E 



(24) 



Here and k n denote the parallel and the perpendicular momenta to plates (in z-direction), respectively. Note 
that, both contributions related to TE mode and TM mode are considered to be the same, hence the final energy 
should be became twice. After the usual Wick rotation and using Eqs. flTTl) and (PMl) with x — x' and carrying out the 
integration over the space then over solid angle in the cZ-dimensional Euclidean space we have 



AR- 



' " 27r)"-ir(^i) J E E ^ 



E 



(25) 



where k E = w 2 + k\ and using the Dirichlet boundary condition on the walls 

k n = , n — 1, 2, 3, ... 
a 

The one-loop photon correction to the vacuum energy in free space is 

AR$ = i [ d 3 x G F (x,x)[-i(g^k 2 - k»k v )6 3 ], 



(26) 



(27) 



where Gf(x,x ) is the propagator of Electromagnetic field in free space in Feynman gauge (£ = 1). After carrying 
out the first integrations and multiplying by f- factor we have 



A< = 



\2aL d - 2 5^^ 

(27r)^r(^i) 



dk F k E 2 



k 2 A- k 2 
dk kE + k 



k% + k 2 + fi 2 ' 



(28) 



Changing of variable as k' = and integrating of kE then expanding these terms about d = 4, all of the divergent 
terms should be eliminated by appropriately adjusting our regulators, for both free and bounded cases. Finally by 
integrating of x we obtain 



AR^-=AR^-AR Fh 



'(i) 



Q{a 2 



3e 2 L 2 /i 2 
4tt 2 



E 



2 , °V 



2 + 2 7 + ln(- 



( 2a \2 
\Lm> 



(29) 



dk'\ k' 2 



a 2 fi 2 



,12 i a 2 At 2 

-2 + 2 7 + ln( f ) 

\Lm) 



+ 0(a 2 ). 



Now, we can use the Abel-Plana Summation Formula (APSF)[5l|, which basically converts a summation into an 
integration, 



E /(») = 



/(0) 



dxf(x) 



dt 



-[f(it)-f(-it)}. 



(30) 



We apply this formula for Eq. (|2l?|) to get (see Appendix for details) 



AE^- 



^e 2 L 2 p? 
47r 3 a 



2 + 2-/ + In 



2 



dt 



Tr 2 t 2 — a 2 /j, 2 



1 



L 2 m 2 /z 2 
4tt 2 

2 + 27 + In 



/ L 2 m 2 



\4a 2 7r : 



/ 2.2 2 2\ 
■(n t — a n ) 



+ 0{a 2 ). 



(31) 



It is obvious that as fi goes to zero, photonic loop does not have any contribution to 0(a) radiative correction to the 
Casimir energy. 
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IV. FERMIONIC LOOP 



In this section, we calculate radiative correction to the Casimir energy due to fermionic loop w ). We use the 

MIT bag boundary condition on the plates. According to MIT bag boundary condition there is no flux of fermions 
through the boundary, this means that 

n„f = 0, (32) 

where j M indicates the current of the Dirac field and is the normal unit vector to the boundary, or more strictly 
it implies to complete confinement of the spinor field. Then MIT bag boundary condition turns out to be (53l - [55| 



(l + i(n.j))iP(x)=0 ) (33) 

which is satisfied on the boundary, more accurate on the two plates. Applying this condition to Dirac spinor field, 
one can derive 

pacot(pa) — —ma, (34) 

which determines quantized modes. Two limits are interesting to calculate, small mass and large mass limits . For 
small mass limit the solutions of Ea. ([34l) are (see for more details [56|) 

p n = (n+-)~ with » = 0,1,2,... (35) 
2 a 

where p n denotes the parallel momenta to the plates (in z-direction) . Now, for the bounded space we have 

/a/2 pa/2 pa/2 

d 3 x(n\n i \n)=i <£ d 3 x + 0(a 2 )=i d 3 x S B (x,x)[i(p6 2 -5 m )}+0(a 2 ), (36) 
-a/2 J— a/2 J -a/2 

where Sb{x,x ) is the propagator of spinor field between plates pSSl l56l] 

S B (x,x) = * fp [ 2 P 2 +m 2 -e-^-'V^-^e-^-A (37) 

a J 2n J (2-K) a ~ l ^ Q w z — p ± — p„ — m z 

Here p± and p n indicate the parallel and the perpendicular momenta to the plates, respectively. After the usual Wick 
rotation and Carrying out the integration, one can obtain 

A P (i)_o r x-2 T ^ 1 />2r, (i^+P^2 mSm 

-* L r(^i) (2tt)^ J Pe 2^pi + pi + m 2 + p E + pl+pi + m i h w 

where p 2 E = w 2 + p\ . Similarly, for the free space we have 

AE F {1) = ijd 3 x S B (x,x)[i(pS 2 - 6 m )]. (39) 

Using Eq.® we obtain 

A£ F (1) = ^"t'' 2 "J, i d P*PE 2 i ~( f E+ f )S \ + a ^ J - (40) 
F (2i) d - 1 T(^) j E J 2np E +p 2 +m 2 p 2 E +p 2 + m 2 ' V ' 
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For the massless case the radiative correction to Casimir energy due to fermionic loop is 



A££ as - = AE^ - AE'( 1] + 0(a 2 ) 



8L 2 n— 



{2n) d ^T{^-) 
,a 2 mfi ( 

+ 2 



dpEPs 2 



1\2 , « Pb \x 
2> 



-i (n -I- 1 ^2 i a Pe i a 2 m- 



dp 



"' 'E 



dp 



P' 2 + 4£ + ^ 



+ 0(a 2 ). 



(41) 



Changing of variables as f/ = |^ and integrating of pe then expanding these terms about d = 4 and /i = 0, as 
mentioned before, all of the divergent terms by appropriately adjusting our regulators, are eliminated. Finally, by 
integrating of x we obtain 



A £Cas. 



e 2 L 2 m 2 



167T 2 a I L ' V 7T^ 

71 = 



E 



I>2 



-6 + 2 7 + In 



V m 3 / 



OO / 9„ 9 



/ a 2 rn 2 , /2 

- 6+27+l,, Hs- 



+ £ 



+ («+ o) 



1\2 



10 - 67 + 61n2 - 3 In 



/ a-rr \2 

U 3 / 



10 - 67 + 6 In 2 - 3 In 



am 1 „/2 



(■ a7T \2 



+ (9(a 2 ). 



Here we need another type of APSF to convert the sum into integral, 

00 1 r°° r°° dt 

Applying this formula in Ea. (|42[) (see Appendix for more details) it becomes 



A£ ;Cas. 



00 Jr2 m 2 



L 2 m 2 y / —a 2 m 2 + ttH 2 



4-4 7 -21n(- 



-a 2 m 2 + n 2 t 2 



8aTr 2 (e 2vt + 1) 

To calculate this expression we use the following summation formula: 



{2air) 2 



) +0(a 2 ). 



OO 1 



(42) 



(43) 



(44) 



The final result of the radiative correction to the Casimir energy due to fermionic loop, for the small mass case, is 
obtained as follows: 



A£; Cas. = y ( l)"e 2 L 2 m 2 / K ^ 2nam) + 2namK 1 (2nam) (j - 2l n7 r - ln(^^) 



16an 2 n 2 

n—L 

The other interesting limit is the large mass limit. In this case, Ea. ([34l) turns out to be 

matan(pa) = —pa. 

Now, the solutions are 

nir 

p n — — with n — 1,2,... 
a 



(45) 



(46) 



(47) 
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The radiative correction to Casimir energy in this case, using Eqs. and (|4T))) . becomes 



AE] 



Cas. 



A£#> - AE'^> + 0(a?) 



'(i) 



8L 2 tt— 



+ 



(27r)"-ir(^i) 

a 2 mS r 



E 



dpEP E 
1 



E 



(n 2 + ^)S 2 



dp 



0(a 2 



(48) 



In a similar way for obtaining Eq. (|42[). now we have 

e 2 i 2 m 2 



A£ ;Cas. 



167r 2 a 



E 



a 2 m 2 i n 2 



a 2 m 2 , 12 



, + n 2 ( - 6 + 2 7 + ln(^- 2 



> +P *( -6 + 27 + m ^ a ; ) 



E 



a 2 m 2 , ^2 

+n 2 (l0-6 7 + 61n2-31n(^-^) 

\ m. ' A t 



, + ^(l0-6 7 + 61n2-31n(^-i- 



0(a 2 



(49) 



Here we can apply the APSF, Eq. (f3"U)) . as we have used in photonic loop, for the above expression, to obtain 

4 - 4 7 - 2 ln( 



A £Cas. = _ 



e 2 L 2 m 3 



32tt 3 



4-4 7 + 21n(^)l - 



dt 



00 e 2 L 2 m 2 v / -a 2 m 2 + 7r 2 t 2 



8an 2 (e M - 1) 



,—a 2 m 2 + TT 2 t 2 



(2air) 



) 



The first term, which appears due to the zero term in APSF, is independent of distance between plates a. Therefore 
this term has no impact on the physics of problem. Ignoring the first term and using the following equation: 



E' 



-2irnt 



1 



we have 



Mi 



Cas. 



E 



(-l)™e 2 L 2 m 2 / . T ^ . ( IQnam 
Ko(2nam) + 2namKi(2nam) ( 7 — 2 In7r — ln( — ) 



16em 2 7r 2 



(51) 



(52) 



Since the photonic correction vanishes up to order a, Eqs. (H5|) and (|5"2"]) are, in fact, our corrections to the well-known 
BRW result for the small mass and large mass limits, respectively. We have compared our results with BRW one 
in Fig. [1] This figure shows that one loop correction due to photonic and fermionic loops is comparable with BRW 
result only when the distance between plates is about one to thirty times of Compton wavelength of fermion. 



V. CONCLUSIONS 



Wc have calculated one-loop radiative correction to the Casimir energy due to photonic and fermionic counterterms 
within the renormalized perturbation theory for QED theory. The topology considered here is two perfectly conducting 
parallel plates in (3+1) dimensions. We have used Dirichlet boundary condition for Electromagnetic field and MIT 
bag boundary condition for electron. To control ultraviolet divergences we have used dimensional regularization and 
a photon mass fj, also is used to control infrared divergences. It is found that photonic loop does not have any 
contribution up to order a. Our result for fermionic loop, up to this order, have had a minor correction to the BRW 
result (341 . 135I ] which they have considered only one diagram of three at this order. We have compared our results 
with the one they have obtained in Fig. [TJ . 
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FIG. 1: The figure shows the O(a) radiative correction to the Casimir energy per unit area due to the fermionic and photonic 
loops vs the distance between plates (in terms of particle Compton wavelength) for two cases: small mass limit (solid line) and 
large mass limit (dashed line). Dot-dashed line is BRW result. Log- plot has been also depicted therein for comparison. 



Appendix: Calculation Of The Branch-Cut Terms 



In this Appendix we calculate two types of branch-cut terms which appear in Eqs. (j!?9")) and (pL2"]) . If we note to 
Ecis. (j29]) and ([42]) . regardless of some coefficients, they are completely the same. Hence we consider only branch-cut 
terms of Eq. 



A£ ,Cas. 



e L 2 m 
I6n 2 a 



E 



i 2 m 2 . 1,„ 

— + ( n + o) 2 



1\2 



-6 + 2 7 + In 



-6 + 2 7 + In 



2a \2 



( 2a \2 



E 



71=0 



1\2 



10 - 67 + 61n2 - 3 In 



^ + (n+i) 



( an \2 



OO / 9„ 9 



10 - 67 + 6 In 2 - 3 In 



a m _,_ „/2 



0(a 2 



For the ease of calculation we use these change of variables: b = —6 + 27 , b' = 10 — 6 7 + 6 In 2, A — ' L 
B = (^) 2 , B' = {f£) 2 and C = We can apply (APSF) 



2 L 2 rn 2 
16ir' 2 a ' 



x -1 />oo />C 

£/(n+-) = / dxf(x)-i \ 

n=0 1 J ° J ° 



dt 



e 2ny + 1 



[f(it)-f(-it)}. 
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Now, the integrand of the last term becomes 

f(it)-f(-it) = Ab(jc+w- v^+f^f) 



A" 



+ Ab' L/C + (it)* - VC + (-H) 2 ) 



,C + (-it) 2 



- aA^c + ^M ^yV ^K-^M " ' y"' )). (53) 

Choosing, C = |C|e* ec , t = \t\e l6t , we have for the first term 

X = Ab\yC + (it) 2 - + (-it) 2 ^ 

'\C\e w c + (e"^|V 2e <) - y/\C\e ie c + ( e - iw |t| 2 e i2fl *)^ 

<«« - e^e^j y/\C\e i ( B c+*- 2B *) + |i| 2 
= 2iAbsm^Vt 2 - C (54) 

where one should note that e J ( 9 c+7r-26i t ) _ _^ This term has been calculated for i > y/C. For i < \/C this term is 
exactly zero. Similarly for second term, we have 

y = VcTWH^P 1 ) vc + HtrH c + i B lt)2 ) 

= - ( \/C + e»f 2 - \JC + e-^t^j In B 

+ (y/C + e"i 2 ln(C + e m t 2 ) -y/C + e-"t 2 ln(C + e^t 2 )^ . (55) 
The first term is exactly X (ignoring the overall factor), so we only need to calculate the second term: 



Y' = \/C + e^t 2 ln(C + e m t 2 ) -\/C + e-™t 2 ln(C + e~ m t 2 ) 

= y/C + e^t 2 \n(e l7r (e- tn C + f 2 )) - y/C + e-^t 2 \n(e~"{Ce" + t 2 )) 

= in(\/C + e™t 2 + VC + e-"t 2 ) + (\/C + e^t 2 - \/C + e -"t 2 ) ln(i 2 - C). (56) 
Now, the first term is again similar to X but with the plus sign between its terms, so we easily obtain 

y/C + e^t 2 + VC + r»i 2 = 2y/t 2 - Ccos |. (57) 



For the second term, using Eg . (|54"]) and Eg . ([ST)) , we have for t < \/C 



y/C + e™t 2 ln(C + e™t 2 ) - y/C + e-^t 2 ln(C + e~ i7r t 2 ) = 2i\/t 2 - C\n(t 2 - C). (58) 
For t < \/C this term is exactly zero. Similarly, third and forth terms in Eq. (|53[) is easily derived. Finally by 
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replacing values of C, B, B', b, b', Eg. (1551) becomes 

e 2 L 2 m 2 \/ —a 2 m 2 + n 2 t 2 



/(it) _ f(- it ) 



8aiT 2 



4 - 47 - 2 In 



-a 2 m 2 + n 2 t 2 



Hence our radiative correction becomes 

,2 r2 m 2 



dt 



e" L 2 m 2 V —a 2 m 2 + ir 2 t 2 
8an 2 (e 2 ^ + 1) 



4-4 7 -2m(- 



(2a7r) 2 
a 2 m 2 + n 2 t 2 



(59) 



(2an) 



) +0(a 2 ). 
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